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Abstract 

While the finite-temperature effective potential in a gauge theory is a gauge-dependent quantity, 
in several instances a first-order phase transition can be triggered by gauge-independent terms. A 
particularly interesting case occurs when the potential barrier separating the broken and symmetric 
vacua of a spontaneously broken symmetry is produced by tree-level terms in the potential. Here, 
we study this scenario in a simple Abelian Higgs model, for which the gauge-invariant potential is 
known, augmented with a singlet real scalar. We analyze the possible symmetry breaking patterns 
in the model, and illustrate in which cases gauge artifacts are expected to manifest themselves most 
severely. We then show that gauge artifacts can be pronounced even in the presence of a relatively 
large, tree-level singlet-Higgs cubic interaction. When the transition is strongly first order, these 
artifacts, while present, are more subtle than in the generic situation. 
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I. INTRODUCTION 



Cosmological phase transitions related to the spontaneous breaking of symmetries in 
fundamental physics are believed to be potentially connected with the deepest questions 
concerning the early evolution of the universe pQ . They may pertain to such diverse topics 
as the origin of seed intergalactic magnetic fields (see e.g. [2]), the excess of baryons over 
anti-baryons (e.g. |3j) and a possible isotropic background of gravitational waves (see e.g. 
[1]). These topics are highly timely, especially in view of the recent claimed detection 
of small but non- vanishing intergalactic magnetic fields [HI E], progress in electric dipole 
moment searches [7J Ej, dark matter [H] and direct collider searches [TO] for signatures of 
electroweak baryogenesis [H] , and, finally, with a new generation of experiments looking for 
gravity waves [12], that will soon boost the already significant results of current detectors 

The possibility of an electroweak phase transition (EWPT) associated with electroweak 
symmetry breaking (EWSB) is especially relevant. In the Standard Model (SM), EWSB 
entails the Higgs field acquiring a non-vanishing vacuum expectation value (vev) that breaks 
the SU(2) xU(l)y gauge group down to U(l) e . m . and generates masses of the weak gauge 
bosons and the SM fermions. The nature of EWSB is governed by the interplay of SM 
gauge interactions and the Higgs quartic self-coupling, which also determines the value of 
the Higgs boson mass, m#. The results of lattice simulations indicate that for m# < 70 — 80 
GeV, EWSB occurs via a first order EWPT, while for a heavier Higgs, the transition is a 
cross-over [Hj. Given the present lower bounds on Mh obtained from LEP, Tevatron, and 
LHC searches [T5] - [T7] . one would conclude that an EWPT would not have occurred in a SM 
universe. On the other hand, extensions of the SM scalar sector can readily lead to a first 
order EWPT as well as associated phenomenology for collider searches. In the context of 
electroweak baryogenesis the EWPT must be strongly first order in order to prevent excessive 
washout of the baryon asymmetry by sphaleron processes. Paradigmatic extensions to the 
Standard Model Higgs sector yielding a strongly first order EWPT include the minimal 
supersymmetric extension to the Standard Model (MSSM) with a light stop [HJ [19], and 
theories (supersymmetric or not) that include one or more extra gauge-singlet fields [201423] 
(other scenarios are also possible, see e.g. [21]). These models typically predict distinctive 
collider signatures in regions of parameter space associated with a strong first order EWPT. 
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Consequently, rapid progress in searches for the SM Higgs at the Large Hadron Collider 
[T7] will soon impact our understanding of a possible EWPT, elucidating whether or not 
the transition was strongly enough first-order for successful electroweak baryogenesis [TO] , 
whether it could have impacted thermal relic densities [251 EE] , and whether it could have 
left any detectable imprint in the diffuse background of gravitational radiation [271 128] . 

Once the field content of a theory is specified, the character of resulting phase transitions 
relies on the computation of an effective potential, V e g, while its dynamics follows from the 
associated effective action, S e g. Although the most robust techniques for computing these 
quantities employ non-perturbative methods, such as discretizing the theory on a lattice 
[TU EHJ [30] , in practice the resulting computational cost makes a perturbative calculation 
by far more feasible and, historically, preferentially pursued for phenomenology. However, 
perturbative calculations of the effective potential generically lead to gauge-dependent re- 
sults, as pointed out long ago by Dolan and Jackiw [31] (see also Ref. [32] for early work 
on gauge-dependence and symmetries at high temperature). Although it is straightforward 
to maintain gauge-invariance when carrying out calculations in the symmetric phase (where 
the vev is gauge-independent), such as computing the temperature of a second-order phase 
transition [33H35] . it is doing so for the broken phase requires additional care. The generic 
dependence of the effective action on the gauge choice is described by the so-called Nielsen 
identities [36] and their generalizations [37]. In practice, the gauge invariance of the effective 
action is guaranteed when the background field <p(x) is an extremal configuration, i.e. one 
that satisfies the equations of motion (in the case of the effective potential, (p(x) = <^ m in 
is the value of the field corresponding to a minimum of the potential) 1 . Typically, gauge 
dependence stems from an inconsistent truncation of the perturbative expansion [38]. This 
leads, in turn, to effects in the critical temperature [39], in the bubble nucleation rate [40J, 
and in the sphaleron transition rate [38j for a first order phase transition, ultimately result- 
ing in unphysical gauge dependence in observable quantities such as the spectrum of gravity 
waves produced by bubble collision or turbulence [28J. 

The problem of gauge dependence as it relates to the description of cosmological phase 
transitions has recently attracted renewed attention. Following earlier studies [JTH33], 

1 More generally, all background fields must be in extremal configurations. For example, the electroweak 
sphaleron involves non- vanishing scalar and gauge fields at the saddle point of the effective action. 
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Ref. |38j addressed the possibility of producing a consistent, order-by-order gauge-invariant 
result in the perturbative expansion of the V e s and S e g. In the SM, this approach yields 
a gauge-invariant critical temperature Tq at one-loop order in a straightforward manner, 
and it can be refined to reproduce leading terms in the "daisy resummation" in a gauge 
invariant manner. Doing so reproduces trends with model parameters that are observed in 
lattice studies, including the dependence of Tq on the top squark soft mass parameters in 
the MSSM. Application to computation of a gauge invariant sphaleron rate is also feasible. 
On the other hand, a gauge-invariant computation of the bubble nucleation rate in the SM 
requires going beyond one-loop order. 

As an alternative to the SM, the Abelian Higgs model provides a theoretically attractive 
"laboratory" in which to assess various approaches to obtaining gauge-invariant quantities 
associated with symmetry-breaking. Apart from calculational ease, this model allows for 
the computation of a gauge-invariant effective potential using a Hamiltonian approach jH] , 
making a comparison with results obtained with other methods possible. In Ref. [28J, we 
explicitly addressed the case of an Abelian Higgs model, and calculated for a full set of 
gauge choices the impact of gauge-dependence on physical observables. Doing so allowed us 
to directly compare the results of the computation for a generic gauge choice with the 
gauge- independent calculation jHJH5]. 

While computationally tractable, the Abelian Higgs model arguably carries limited phe- 
nomeno logical interest. Among the more relevant SM extensions mentioned above - such 
as MSSM with a light stop or a gauge-singlet extension to the Higgs sector - it is possible 
to make the electroweak phase transition strongly first order via interactions that are gauge 
independent (though the full Kff remains gauge-dependent). The simplest cases involve 
introduction of tree-levels term in the potential of the type SH^H or S 2 H^H in the case 
of the real singlet extension. The tree-level cubic operator can produce a large potential 
barrier between the broken and unbroken phase at the electroweak phase transition, while 
the quartic interaction may allow for a lowering of Tq in a manner compatible with collider 
constraints on ma- As these operators are manifestly gauge- invariant, one may inquire as 
to whether perturbative computations of the EWPT properties in the associated models are 
quantitatively less susceptible to gauge-dependent artifacts than in either the SM or in the 
Abelian Higgs model. Indeed, Refs. [TU1 HH] have recently suggested that such a situation 
may occur. 
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In what follows, we study the issues described above in some detail. For the sake of 
comparing with a known, simple gauge-independent result, we shall again use the Abelian 
Higgs model, supplemented with a gauge-singlet real scalar field (which does not impact the 
gauge-dependence structure of the theory) and retaining only the tree-level cubic operator, 
SH^H . Arguably, this model is the simplest prototypical electroweak-like theory that can 
exhibit a strongly first order phase transition driven by tree-level cubic terms. In view 
of the recent LHC results pointing to a relatively heavy Higgs mass, singlet extensions to 
the electroweak scalar sector have additionally become phenomenologically more appealing, 
making an assessment of the gauge artifacts in the effective potential even more timely. We 
show that such gauge artifacts may arise even in the presence of a large tree-level singlet- 
Higgs cubic coupling. However, we also find that the gauge-dependence is less pronounced 
when the tree-level and loop-induced cubic interactions conspire to generate a sizeable barrier 
between the broken and unbroken phases at low temperatures. 

The remainder of this paper is organized as follows: the next section describes in detail the 
theory we study, including explicit calculations of the gauge dependent terms in the effective 



potential. The following Section HI gives an outline of the possible patterns of spontaneous 
symmetry breaking, and describes the effects of gauge choices on various quantities of interest 
(including the critical temperature, the latent heat, and a measure of the strength of the 



phase transition). Finally, Section IV summarizes and concludes. 



II. THE ABELIAN HIGGS MODEL PLUS A SINGLET SCALAR 

We examine the gauge dependence of a simple Abelian Higgs model containing a single 
complex scalar $ charged under a local U(l) gauge group, and a real scalar singlet field s. 
For our purposes here, we consider only a cubic coupling between the two fields, as the latter 
can generate a tree-level barrier between the broken and unbroken phases and, thus, can 
increase the strength of the phase transition. The effect of the quartic operator discussed 
above is more subtle, so for simplicity we focus on the cubic interaction. At tree-level, the 
potential is 

V {*, s) = IX^®) 2 + + |A 2 s 4 + \£s 2 + lEs&®. (1) 
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It is useful to separate $ into real and imaginary parts ($ = h + ih!) and then rotate into a 
basis such that only the real part gets a vev. 

Rather than specifying the 5 coefficients explicitly, we find it more convenient to specify 
the vevs, the tree-level mass eigenstates, and the mixing of the mass eigenstates, and use 
this to set the tree-level potential. The tree-level mass-squared matrix is 

Uk> + A + Es Eh \ 
" I Eh 3\ 2 s 2 + fi 2 J 

/ens/? — sin 9 \ I m? \ / cos/? sin 9 \ 

(3) 




cos 9 sin 9 
— sin 9 cos 9 



where the rotation angle 9 gives the mass eigenstate mixing: 



2M- 



tm(2 «) = Mf^h- < 4 > 

The quantities mi, rri2, and 9 are most relevant to collider phenomenology, as they determine 
production cross sections and decay branching ratios. In particular, a non-zero mixing angle 
9 can weaken the collider constraints on the lightest mass eigenstate since its effective cou- 
pling to gauge bosons is reduced by cos 9. This effect opens up the possibility of a (stronger) 
first order phase transition by allowing for a smaller Higgs quartic self-coupling, although 
we do not study this effect in detail here. The corresponding minimization conditions are 
given by 

rfV 

' (Ai/i 2 +nl + Es)h = 0, (5) 



dh 

^ = A 2 s 3 + fis + \Eh 2 = 0. (6) 



From these we solve for Ai, A 2 , ^1, an d E in terms of vacuum expectation values of h and 
s, masses mi and m 2 , and the angle 9, assuming that both h and s are non-zero. Note that 
for E < (so that the (s) > 0) and m\ < m 2 , we require < 9 < 90°. Also, A 2 is negative 
for sufficiently large E, so not all values of 9 lead to stable potentials. 

For h = 0, extrema occur at s = and — for fi 2 , < — at s = ±|/i 2 |/\/A 2 . An additional 
three extrema can occur for h > 0, whose locations are trivially determined by solving the 
cubic equation in s obtained from from combining Eqs. ^ and ([6]): 

A 2 s 3 + {A ~ ~ ^ = °- (7) 



FIG. 1: Contours of the tree-level potential for mxjmi = 0.4 and six different values of 9. Red 
(blue) contour lines denote higher (lower) values of the potential. The Higgs and singlet fields vary 
along the horizontal and vertical axes, respectively. The origin, which is in the center of each plot, 
is a maximum in (a), (b), and (f), a saddle point in (c) and (e), and a minimum in (d). 

Since Eq. ^ has no term proportional to s 2 , at most two of the extrema can be in any one 
quadrant. Also, since there are no linear or cubic terms in h, all maxima must lie along the 
s-axis. Using this knowledge, one can enumerate all of the different combinations of minima, 
maxima, and saddle points to obtain all of the different possible potential types. 

we observe that the location and character of the extrema depend on 
four independent parameters. For example, by scaling out a factor of Ai we may take these 
parameters to be [if/Xi, f^l/Xi, A 2 /Ai, and E/\\. We may trade two of these parameters for 
one each of the non-zero vacuum expectation values of s and h, respectively. The remaining 
two parameters then determine 9 and the ratio of masses m\jm,2- The depth of the potential 
at one of the minima is then fixed by the fifth remaining parameter in the potential, which 
we can trade off for one of the masses. Therefore, we can keep two vevs and one of the masses 
fixed and just vary 9 and the ratio m 1 /m2 to explore all potentials that are not related by an 



From Eqs. (5|6 



overall rescaling. Fig. [T] shows six different representative potentials with constant 777.1/777.2. 

Panel (a) shows the potential with a small positive value of 9, corresponding to a small 
negative E. There is a minimum in each quadrant of the h-s plane separated by four saddle 
points with a maximum at the origin. Increasing the angle 9 (panels (b) and (c)) merges 
some of these features onto the s-axis. First, the two minima at s < merge, then the two 
saddle points near s = merge onto the origin such that the origin is no longer a maximum. 
Increasing 9 further (panel (d)), pushes the minimum along the s-axis up to the origin. At 
this point there is a tree-level barrier between the broken and symmetric phases. In panel 
(e), this barrier disappears and there are only the two electroweak minima and a saddle 
point at the origin. Finally, panel (/) has a small value of \E\, but large enough to destroy 
the two metastable minima in panel (a). A further rotation to 9 = 90° would reproduce a 
potential of a type similar to what shown in panel (a). This basically exhausts all of the 
possibilities with symmetry-breaking minima: the only other potential type occurs at larger 
ratios of 7712/777,1. It is similar to type (d) except that the h = minimum splits into two 
minima and a saddle point along the s-axis. 

A. Quantum corrections 

The one-loop corrections to the effective potential are given by 

where n» is the d.o.f. for each particle, m 8 is the particle mass, q is the renormalization 
scale (which we set to 1 TeV), and c = 1/2 for transverse gauge boson polarizations and 
3/2 for all other particles. The Higgs and scalar masses are given by the eigenvalues of the 
tree-level mass matrix (Eq. [2]). The three gauge boson polarizations each contribute a mass 
m 2 auge = g 2 h 2 . We focus on gauge 2 , in which there are gauge-dependent masses for the 
Goldstone boson and the ghost: rn 2 old = Ai/t 2 + /j,1 + Es + £g 2 h 2 and m 2 ghost = £g 2 h 2 , with 
nghost = —2. There is an additional degree of freedom from the gauge boson's unphysical 
time-like polarization which exactly cancels one ghost degree of freedom. In all of our 

2 In a previous study , we compared gauge to the gauge- independent potential in Ref. [H] which uses 
a Hamiltonian formalism. The latter is more computationally intensive, and it tends to closely resemble 
Landau gauge (£ = 0) in the Abelian Higgs model. 
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models the first derivatives (dVi/dh and dVi/ds) are negative at the tree- level minimum. 
This pushes the vevs further away from the origin than they were at tree-level. 
At finite temperature, the one-loop corrections are 



2tt 

where 



v -^rwW^ 



J^2 



(9) 



J(x 2 ) = dy y 2 log (l - e" v^+^j . (10) 



In the high-temperature (low-x) limit, 



J( x 2) ~ + IL^ - -X 3 - — tog - - 0( 6) ( , , ) 

v ; 45 12 6 32 a 6 v ; v ; 

where loga& = | — 27^ + 21og(47r) and 7^ is the Euler constant |Hj. J(x 2 ) is not analytic 
at x 2 = 0, and it is complex for x 2 < 0. In our calculations, we approximate J(x 2 ) with a 
cubic spline, taking only the real component for x 2 < 0. 

At high temperature, the validity of the perturbative expansion of the effective potential 
breaks down. Quadratically divergent contributions from non-zero Matsubara modes must 
be re-summed through inclusion of thermal masses in the one-loop propagators [48, 49J: 
m 2 (<j)) — > m 2 e ff((f)) = m 2 (<f)) + mf fterm (T). This amounts to adding thermal masses to the 
scalars and gauge boson longitudinal polarizations: 



Ai/3 + 2 / 4 
A 2 /4 



M 2 ^M 2 +T 2 ( / ] (12) 



m 2 gold ^ m 2 gold + T 2 (X 1 /3 + g 2 /4) (13) 

—gauge 

Note that the coefficients of the T 2 terms are ^-independent. 

When the phase transition is second-order or very weakly first-order, or when the temper- 
ature is very high, even the re-summed potential may not be reliable. Loops that are either 
infrared divergent or dominated by the infrared regime contribute linearly in temperature 
and can ruin the perturbative expansion (see, e.g., Ref [50J). Each additional such loop 
contributes roughly XT/M, where A is the relevant coupling and M is the relevant mass 
scale. Substituting the gauge boson mass for M and g 2 for A, we see that the perturbative 
expansion should hold as long as h/T > g. We warn the Reader that for certain parameters 
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in what follows this criterion breaks down. As a result, the one- loop expansion might have 
limited validity in those cases. This is particularly true for cases 1 and 2 with g = 0.5 
and case 3 with 6 > 70°. However, it is important to note that one could lower the phase 
transition temperature by e.g. extending the gauge group and adding extra gauge bosons. 
Extra degrees of freedom enhance the finite-temperature contributions relative to the tree- 
level potential, so symmetry breaking happens at lower temperatures. This in turn would 
make the one-loop perturbative expansion more reliable, without qualitatively changing the 
nature of the explicit gauge dependence. Since the appearance of the one-loop gauge de- 
pendence is not tied directly to the perturbative validity, and since our primary interest is 
in providing proof of existence for the gauge-dependence issues, we leave the perturbative 
breakdown problem to future studies. 

III. PATTERNS OF SPONTANEOUS SYMMETRY BREAKING 

We now study the gauge dependence in a few representative models with different quali- 
tative features. In particular, we examine the gauge dependence of the nucleation tempera- 
ture T* and various measures of the phase transition strength. Nucleation occurs when the 
three-dimensional action S3 of a nucleated bubble satisfies pa 140 (see Refs. [SI1E2] for 

original work on cosmological phase transitions), and we use this criterion to define T*. The 
phase transition strength has often been characterized by A(f>, the jump between the vevs of 
the two phases. However, as noted above, this quantity is ^-dependent. Alternate, physically 
meaningful measures include (a) a, the difference in energy densities at the two vevs , and (b) 
a measure of the phase transition duration defined as /3/H* = T*(d/ 'dT)(S%/T) where 

is the Hubble parameter at the time of the transition. When the phases are degenerate 
and there is no supercooling, a is equivalent to the transition's latent heat. The quantity f3~ l 
effectively measures the strength of the transition, with f3 = 00 for a second-order transition. 

We use the Cosmo Transitions package [53] to determine the phase structures and calculate 
the nucleation rates. 
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TABLE I: Model parameters for three illustrative cases. 





{htree) (stree) 


mi 777.2 


e 


E 


Case 1 


240 GeV 48 GeV 


120 GeV 35 GeV 


-i° 


-0.96 GeV 


Case 2 


240 GeV 48 GeV 


120 GeV 180 GeV 


45° 


-37.5 GeV 


Case 3 


174 GeV 174 GeV 


30 GeV 75 GeV 


60° 


-11.8 GeV 



A. Gauge fields driving transitions 

First, we examine two cases in which the transition would be second-order without the 
inclusion of massive gauge fields. Table [I] contains the corresponding model parameters. The 
two cases are quite similar: they have the same vevs, and they both have a saddle point at 
h = s = with no other extrema along the s-axis (see Fig. [ije)). However, case 1 has a 
very small cubic term while in the second case E is of the same order as the mass scale of 
the theory. The gauge dependence exhibited in case 1 is therefore unsurprising: since it has 
relatively weak coupling between the Higgs and the scalar singlet, we expect it to show the 
same sort of gauge dependence as an uncoupled Abelian Higgs model [28]. This is indeed 
the case, as seen in Fig. [2] (green lines). 

For low values of the gauge coupling g, the gauge dependence is quite pronounced. The 
initial symmetry breaking (solid lines for the two cases) is weakly first-order (e.g., a/T^ <C 1) 
from £ = to £ « 1.5. At higher £, a second-order transition initially breaks the symmetry 
(dashed lines in Fig. [2]). The first-order transition then proceeds from a high-temperature 
broken phase to the low-temperature broken phase at larger values of h and s. Above 
£ ~ 3, the barrier between the two broken phases disappears and there is no first-order 
phase transition at all. 

For higher values of g, the gauge dependence is not as pronounced. The phase transition 
is much more strongly first-order, with a/Tf roughly a factor of 10 higher than it is for 
low g. The initial symmetry-breaking still turns second-order at high £, but the subsequent 
first-order transition persists up to £ = 5 with about a factor of 2 drop in a. 

In all cases, the appearance of a strong first-order phase transition is associated with 
a large ratio of field-dependent heavy degrees of freedom (in this case, gauge bosons) to 
the Higgs mass. A small Higgs mass decreases the depth of potential (that is, Vq(0) — 
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320 
300 
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260 
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Case 2 — — 


1 _ 
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1 1 1 


1 ~ 





~5" 



10° 



lo- 



ur 



FIG. 2: Gauge dependence in cases 1 (thin green lines) and 2 (black lines) as a function of the gauge 
parameter £. The left panels have a small gauge coupling g = 0.5, while the right have g = 1.0. 
Dashed lines represent second-order symmetry breaking transitions, which may be followed by a 
first-order transition at lower temperature. The marks along the right side of each panel show 
the corresponding quantity calculated using the gauge-invariant method of Ref. |38j . The thicker 
marks include a gauge-invariant treatment of the thermal masses; the thin marks ignore them. 

Vo(v) decreases), while heavy additional field-dependent masses yield larger contributions 
to the thermal effective potential (increasing Vi(v,T) — Vi(0, T)). Both effects decrease 
the critical temperature and increase A0. The presence of the additional field-dependent 
masses decreases the critical temperature because, for a given value of T, V\(v, T) — Vi(0, T) 
is larger for larger gauge couplings and the two minima are degenerate at lower temperatures. 
It increases the value of A0 because dJ/dx — >■ + as x — > oo, so when x = m/T = gh/T 
is large, dV\/dh is small and the vev does not decrease much from its tree-level value. An 
increase in A0 tends to increase both a (a larger separation between phases implies a larger 
difference in mass spectrums, entropy, and therefore latent heat) and (since S3 scales as 
(A0) 3 ). One can achieve a strongly first-order phase transition even for a heavy Higgs, as 
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long as it is somewhat light compared to the other field-dependent masses. 

Interestingly, case 2 (thick black lines) shows almost exactly the same gauge dependence 
as case 1, even though it has a non-trivial cubic term. The important point is that, although 
large, the cubic term is not large enough to cause a first-order phase transition without 
additional bosons that have large couplings to the Higgs field. 

We compare the gauge-dependent calculations to explicitly gauge-independent calcula- 
tions (denoted by marks along the right side of each panel in Fig. [2|. At one-loop order and 
without the added thermal masses, the gauge-independent calculation is simply the value 
of the potential evaluated at the tree- level minimum <&o, where ghost and goldstone degrees 
of freedom exactly cancel. Thermal masses spoil the cancellation, but one can still obtain a 
gauge-invariant result by evaluating the cubic terms in the ring-improved effective potential 
at the tree-level high-temperature minimum, where the tree-level high-temperature potential 
is the same as Vo($) but with thermal masses added to fif and This is the lowest-order 
approximation used by Ref . j3H] ■ Since the potential is evaluated at two different minima, a 
gauge-invariant ring-improved A0 is not well-defined and is not plotted in Fig. [2] 

Ignoring thermal masses, one can see that the gauge-invariant critical temperature must 
be lower than the gauge-dependent critical temperature for any value of £: the gauge- 
dependent critical temperature is defined as the temperature at which V($ m i n , T) = V(3> = 
0,T), but since <& m in is the minimum of the potential, V(§ ,T) > V(Q min ,T) = V($ = 0, T) 
and the gauge-invariant critical temperature must be lower. Conversely, the latent heat tends 
to be larger in the gauge-invariant method. The energy density decreases with increasing 
particle masses, so, as long as the masses are larger at <3>o than at $ m j„ (which is the case for 
weakly first-order transitions), the difference in energy densities between the symmetric and 
broken phases will be larger when evaluated at $o than at Q m in- The addition of thermal 
masses tends to enhance both of these effects. 

The gauge-invariant method produces quite different results from the gauge- dependent 
calculation when the latter predicts a very weakly first-order transition. This is to be 
expected: the two methods perform calculations at very different field values when the 
gauge-dependent A</> is small. When g — 1 and the transition is more strongly first-order, 
the two methods agree much more closely. However, including thermal masses worsens the 
agreement. When m gauge /T > 1, higher-order terms in the effective potential dominate and 
using only the cubic term for ring-improvement is unreliable. The ring-improved gauge- 
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\E\ (GeV) 

5 10 15 20 25 30 35 37.5 
340 I i i i 1 1 1 1 1 




5 10 15 20 25 30 35 40 45 
9 (degrees) 



FIG. 3: Gauge dependence in case 2 with g = 0.5, but with 9 (or equivalently E) varied and £ 
held fixed. The transition is second-order for £ > 3 for all values of 9. 

invariant calculation should not be trusted in this case. 

Fig. [3] explicitly shows the (non-)importance of the cubic term (via 9) in these scenarios. 
Regardless of whether the cubic term is large or small, the basic pattern of gauge-dependence 
is about the same. The phase transition grows more weakly first-order (a decreases) for 
increasing £ for all values of 9. At £ > 3, there is only a second-order transition. Note that 
the transition is most strongly first-order when tan 9 = (h tr ee) / (stree), but still second-order 
for £ > 3. 



B. Tree-level terms driving transitions 

Here we examine a scenario in which the cubic term is critical in determining the strength 
of the phase transition. Superficially, case 3 appears similar to case 2. Both have relatively 
large cubic terms, and both have the topology shown in Fig. However, a small change in 
model parameters can turn the saddle point in case 3 into a tree-level minimum (Fig. [l]((i)), 
thus creating a potential barrier at zero temperature for which the tunneling rate may never 
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be large enough to penetrate. Even without a tree-level barrier, the cubic term is sufficiently 
prominent to create a barrier at relatively low temperature: there is no barrier at T = 0, 
but there is a barrier by T « 100 GeV for 9 = 60°. Slightly smaller values of 9 decrease this 
temperature drastically. The crucial distinction between cases 2 and 3 is that, even though 
both have large cubic terms, only in case 3 is the lowest eigevalue of the mass-squared matrix 
both negative and sufficiently small in magnitude (i.e., a small negative value of /if) that 
the origin can become a minimum at low temperature with the cubic term providing the 
separation between the symmetric and broken phases. 

Since the phase transition is strongly first-order even for g = (that is, without any 
gauge bosons at all), the gauge dependence is not nearly as severe as in cases 1 and 2. It 
is still present though. For example, at 9 = 60°, T* increases by 1.6% from £ = to £ = 5, 
and a increases by only 0.1%. 

By increasing 9 from 55° up to 90°, one moves successively through topologies (d), (e), 
(/), and (a) in Fig. [TJ At around 9 = 70° and \E\ = 9 GeV, the cubic term is small enough 
so that a first-order phase transition requires g > 0. At this point, the gauge dependence 
becomes much more obvious, as is seen in Fig. |4} At high enough 9 and low enough E, the 
symmetry-breaking transition is second-order for all plotted values of £. 

IV. DISCUSSION AND CONCLUSIONS 

The main conclusion of the present study is that the inclusion of a singlet scalar degree 
of freedom does not generally alleviate the gauge-dependence problem in the electroweak 
phase transition, even when it has large couplings to the Higgs. Moreover, a significant, 
tree-level cubic singlet-Higgs interaction does not in itself guarantee a strongly first-order 
phase transition. On the other hand, when the phase transition is strongly first order, the 
gauge-dependence appears to be less pronounced than in the generic case. Such a situation 
occurs either when the gauge coupling is relatively large or when the tree-level singlet-Higgs 
cubic term acts in concert with small negative mass-squared values at h = s = to create a 
potential barrier at low temperature. Otherwise, when the phase transition is only weakly 
first-order, or borderline strongly first-order, the gauge- dependence can be drastic regardless 
of the presence of a cubic term. This dependence may change not only the strength of the 
phase transition, but also its overall character. In such circumstances, one cannot make 
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FIG. 4: Similar to Fig. [3J but for case 3 with g = 0.3. At 9 > 70°, the symmetry- breaking 
transition is second-order for £ = 3 and 5 (dashed lines), followed by a weakly first-order transition. 
At 9 > 80°, the symmetry-breaking transition is second order for all plotted values of £ and the 
first-order transitions have disappeared. The thick/thin black dotted lines show the gauge-invariant 
critical temperature and latent heat calculations with/without thermal mass corrections. 

a gauge-independent determination of whether the transition is first or second-order, nor 
can one even determine whether or not the transition comes from a symmetry-preserving 
vacuum. The explicitly gauge-independent calculation of the critical temperature Tc using 
the lowest-order result in Ref. [38J can give a rough estimate of the transition temperature 
T*, but only when the amount of super-cooling is small 3 , which is hard to achieve with a 

3 In this case, the onset of nucleation occurs for T very close to Tq 
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tree-level barrier. A similar gauge- independent calculation of a is only reasonable when the 
transition is already known to be strongly first-order. 

As emphasized in our earlier work [28J, the appearance of gauge-dependence in physical 
quantities such as T*, a, and should engender caution when attempting to draw phe- 
nomenological conclusions from computations performed in a specific gauge. In the ideal 
situation, a gauge-invariant computation using non-perturbative methods would be used to 
explore various Standard Model extensions that may lead to a first order electroweak phase 
transition, though a comprehensive exploration is at present prohibitively expensive. In the 
meantime, various gauge-invariant perturbative techniques, such as the loop expansion [38] 
or Hamiltonian formulation [33], may at least point to regions of parameter space in a given 
model where transitions of different character occur. If, as we find for the Abelian Higgs plus 
singlet model (and as, perhaps, maybe speculatively be a more general pattern) the gauge- 
dependence of conventional perturbative computations is mitigated by a strongly first order 
transition triggered by gauge-independent terms, one might expect to find rough agreement 
with the results of manifestly gauge-invariant analyses. The present study shows, however, 
that such a conclusion should be carefully qualified on a case-by-case basis. 
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